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SUMMARY 

This research program investigates boundary-layer receptivity in the leading- 
edge region for bodies with blunt leading edges. Receptivity theory provides the link 
between the unsteady disturbance environment in the free stream and the initial amp- 
litudes of the instability waves in the boundary layer. This is a critical problem which 
must be addressed in order to develop more accurate prediction methods for 
boundary- layer transition. The first phase of this project examines the effects of 
leading-edge bluntness and aerodynamic loading for low Mach number flows. In the 
second phase of the project, the investigation will be extended to supersonic Mach 
numbers. Singular perturbation techniques are utilized to develop an asymptotic 
theory for high Reynolds numbers. In a parallel, closely coordinated effort, the 
author is collaborating with Dr. T.B. Gatski of NASA LaRC in utilizing his Navier- 
Stokes code for computations of leading-edge receptivity. 

In this first year of the project, we have developed the asymptotic theory for 
leading-edge receptivity in low Mach number flows. The case of a parabolic nose is 
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considered, with 5 = u)r n /U <X) = 0(1). The asymptotic theory for the unsteady flow 
utilizes the small parameter e = . In the outer region away from the body 

surface, an inviscid solution holds. The viscous analysis in the boundary layer 
involves two streamwise regions. In the region wx/U^ = 0(1), the disturbances satisfy 
the linearized, unsteady boundary- layer equation (LUBLE). Farther downstream 
where uix/U <*, = 0(f J ), an asymptotic form of the Orr-Sommerfeld equation (OSE) 
applies. The inviscid solution in the outer region has been calculated, and the asymp- 
totic eigenfunctions for the LUBLE have been derived and matched with the 
Tollmien-Schlichting (TS) wave of the OSE. This matching shows that the coefficient 
C l of the first asymptotic eigenfunction of the LUBLE determines the amplitude of 
the TS wave. To complete the asymptotic analysis, the coefficient C : of the asymp- 
totic eigenfunction must be determined from a numerical solution of the LUBLE. We 
are presently addressing this task. 

Substantial progress has also been made on the Navier-Stokes computations. 
Analytical solutions for the steady and unsteady potential flow fields have been incor- 
porated into Dr. Gatski's code, greatly expanding the types of free-stream distur- 
bances that can be considered while also significantly reducing the computational 
requirements. The time-stepping algorithm has been modified so that the potential 
flow perturbations induced by the unsteady pressure field are directly introduced 
throughout the computational domain, avoiding an artificial "numerical diffusion" of 
these from the outer boundary. In addition, the start-up process has been modified 
by introducing the transient Stokes wave solution into the downstream boundary con- 
dition. These modifications have significantly improved the performance of the code, 
and we anticipate that within a few months computational results will be available for 
comparison with the asymptotic theory. This will be the first comparison of asymp- 
totic and Navier-Stokes results for leading-edge receptivity. 


- 3 - 


PROGRESS REPORT 

The accurate prediction of boundary- layer transition and the need to control 
this phenomenon are becoming increasingly important in the development of 
advanced aeronautical vehicles. These requirements have focused renewed attention 
on the physics of the transition process. It was known for many years that boundary- 
layer transition is significantly influenced by the free-stream disturbance environ- 
ment. However, the physical mechanisms by which energy is transferred from the 
long-wavelength, free-stream disturbances to the short-wavelength, instability waves 
were not understood. This came to be known as the receptivity problem. 

Recent work by Goldstein and extensions by the author and his students have 
shown that the receptivity process takes place in regions of the boundary layer where 
the mean flow exhibits rapid changes in the streamwise direction. This occurs (a) 
near the body leading edge and ( b ) in any region farther downstream where some 
local feature forces the boundary layer to adjust on a short streamwise length scale. 
The rapid streamwise adjustment in these regions requires that nonparallel mean flow 
effects be included at leading order, in contrast to the parallel flow or slowly diverg- 
ing flow assumptions utilized in stability theory. The current status of research in 
receptivity is summarized in several articles contained in Ref. 1. 

Previous theoretical research in receptivity near leading edges has focused on 
the Blasius boundary layer on a semi-infinite, zero-thickness plate. The motion has 
been assumed to be incompressible and two-dimensional. Goldstein 2 utilized a high 
Reynolds number asymptotic approach to develop a theoretical framework for leading 
edge receptivity. The small parameter « utilized in the analysis is defined by 

2 iJq 

e = {uiv/Ug^ ; note that e 6 is the frequency parameter F for classical stability theory, 
which can also be interpreted as the inverse of a Reynolds number based on the con- 
vective wavelength U 0 0 /w of the unsteady motion. 


Two distinct streamwise regions arise in the asymptotic analysis. Near the 
leading edge where wx/U^ =0(1), the unsteady motion is governed by the linearized 
unsteady boundary layer equation (LUBLE). In the second streamwise region farther 
downstream where = 0(e~ 2 ), the unsteady motion is governed by the classical 

large Reynolds number, small wavenumber approximation to the Orr-Sommerfeld 
equation (OSE). The Tollmien-Schlichting (TS) wave in the latter region corresponds 
to the triple deck structure of Smith 3 . Goldstein examined the asymptotic matching 
of these two regions, and showed that Lam and Rott's first asymptotic eigenfunction 4 
for the LUBLE matches onto the TS wave of the OSE. Thus, the first asymptotic 
eigensolution of Lam and Rott is the precursor of the TS wave, and the coefficient C l 
of this asymptotic eigensolution determines the amplitude of the TS wave. Hence, we 
call C x the "Receptivity Coefficient." 

The available results for the Receptivity Coefficient are restricted to the low 
Mach number, Blasius boundary layer on a semi-infinite, zero-thickness plate. The 
first calculation of a Receptivity Coefficient was presented by Goldstein, Sockol and 
Sanz 5 for the case of an acoustic wave propagating downstream parallel to the plate 
surface. The author and Heinrich 6 ' 7 utilized information on the analytic structure of 
the solution to develop a much more accurate approach for calculation of the Recep- 
tivity Coefficient. Using this improved method, they calculated Receptivity Coeffi- 
cients for a wide range of free-stream disturbances including obliquely incident 
acoustic waves, convected gusts (the linear representation of free-stream turbulence) 
of various orientations, and a von Karman vortex street passing above the plate 
surface. The Receptivity Coefficients for obliquely incident acoustic waves were 
found to be an order of magnitude larger than those for other free-stream distur- 
bances. These large receptivity levels arise arise due to the strong diffraction field 
produced by the interaction of an oblique acoustic wave with the plate leading edge. 
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In order to obtain predictions of relevance to practical applications, it is essen- 
tial that leading-edge bluntness be incorporated into the theory for leading edge 
receptivity. All bodies of practical interest have a finite leading-edge radius r„, 
which enters the receptivity problem through the nondimensional parameter 5 = 
wr n/U oo- Th^ parameter is typically 0(1) or larger in aerodynamic applications, and 
it may even be large enough in laboratory experiments to cause significant deviations 
from the zero-thickness result. In the Receptivity Panel at the ICASE/LaRC 
Workshop on Instability and Transition (see Ref. 1), both the author and Morkovin 
identified leading-edge thickness as an important area for future research. 

A second disadvantage of the zero- thickness geometry is that the influence of 
mean aerodynamic loading cannot be investigated, since any asymmetrical mean flow 
component in the vicinity of the leading edge leads to boundary layer separation. It is 
well known that aerodynamic loading significantly influences both the mean boundary 
layer development and the stability of the boundary layer. Hence, it is reasonable to 
anticipate that aerodynamic loading is an important parameter in the receptivity 
problem as well. In fact, in his laboratory experiments Leehey 8 found that small 
changes in the "angle of attack" of his plate produced dramatic changes in the recep- 
tivity to acoustic waves. Thus, in order to obtain results which are relevant to practi- 
cal applications, it is essential to incorporate both leading-edge bluntness and aerody- 
namic loading in the receptivity theory. 

The present research program extends the asymptotic theory for leading edge 
receptivity to incorporate the effects of leading-edge bluntness, aerodynamic loading 
and compressibility. In the first phase of the program, the effects of leading-edge 
bluntness and aerodynamic loading are being investigated for the low Mach number 
case. In the second phase of the program, the theory will be extended to high Mach 
number flows. The low speed problem is viewed as an essential stepping stone in the 
development of the theory, and also has relevance to the high Mach number case 



where a strong shock produces low speed flow near the nose of the body. 

As an integral part of this research program, we are also collaborating with Dr. 
T.B. Gatski of NASA LaRC in Navier-Stokes computations for leading-edge recep- 
tivity. These computations are an important complement to the asymptotic theory, 
which provides valuable information regarding the dominant physical mechanisms and 
relevant scaling parameters, but may not provide sufficient quantitative accuracy for 
parameter ranges of practical interest. In particular, the Navier-Stokes computations 
allow finite Reynolds number effects to be assessed. The computations can also be 
utilized to investigate nonlinear effects due to free-stream disturbances of finite amp- 
litude. 

For the low Mach number study which forms the first phase of this research 
project, we consider conventional airfoil shapes which have a parabolic leading edge. 
The streamwise length scale for the leading-edge receptivity analysis is the convected 
wavelength U^/w, which is generally small compared to the airfoil chord b. Hence, it 
is only the airfoil characteristics near the leading edge which enter the analysis, and 
the geometry can be restricted to a semi-infinite parabola with the nose radius r n of 
the airfoil. For a symmetric airfoil at zero angle of attack, the flow past the leading 
edge is symmetric. However, aerodynamic loading due to camber or a nonzero angle 
of attack produces a flow asymmetry. In the local leading-edge region, this asymme- 
try corresponds to a flow component around the leading edge from the pressure side 
to the suction side. The application of the method of Matched Asymptotic Expan- 
sions 9 for a nose radius r n small compared to the airfoil chord b shows that, in the 
local leading-edge region, the complex potential for the inviscid mean flow has the 
asymptotic expansion 


+ VooBeffb 1 ^ [Vr - /'A, /2] + ... (1) 


W(z) - $ + ft - [Vr - /V^T/2] 2 
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where z = x + iy. The first term corresponds to the symmetric component of the local 
flow past the leading edge, while the second term is the asymmetrical component 
produced by the airfoil camber and angle of attack. An expression for a e ^ in terms 
of the airfoil angle of attack and camber distribution is given in Ref. 10. In our work 
to date, we have concentrated on the case of leading-edge bluntness for a symmetric 
mean flow (a e ^ = 0). The extension to include nonzero values of a e ^ is conceptually 
straightforward and will be carried out later. Details of our progress on the Navier- 
Stokes computations and on the asymptotic theory are discussed separately below. 

Navier-Stokes Computations 

In the first few months of the project, we focused on the collaboration with Dr. 
Gatski in Navier-Stokes computations for leading-edge receptivity. Dr. Gatski's 
code 11 is well suited to the study of leading-edge receptivity, since it solves the 
Navier-Stokes equations in parabolic coordinates. Hence, no discontinuities in the 
geometry, grid generation functions, or related derivatives are present in the numeri- 
cal approach. Receptivity is a sensitive phenomenon and such discontinuities could 
easily lead to spurious "numerical sources" of receptivity. 

The Navier-Stokes computation of the unsteady flow past a parabola is a diffi- 
cult problem on which previous investigators have had only partial success. Murdock 
attempted an incompressible, Navier-Stokes calculation of receptivity at the leading 
edge of a parabola, and presented a preliminary report in a conference paper 12 . 
However, Murdock never published this work in archival form because his down- 
stream boundary condition caused difficulties which he never satisfactorily resolved 
(private communication, June 1990). Reed et al 13 have presented incompressible 

computations for geometries consisting of an elliptical nose attached to a flat plate. 

c 

Unfortunately, her results appear to be dominated by localized receptivity 14 at the 

t 
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junction between the elliptical nose and the flat plate, and may also be contaminated 
by "spurious sources" related to the numerical treatment of the abrupt discontinuity in 
curvature at the junction between the elliptical nose and the flat plate. 

Both Murdock and Reed introduced the flow unsteadiness by specifying a 
uniform axial pulsation on the outer boundary of the computational domain far from 
the body surface. This boundary condition is an attempt to model receptivity to a 
parallel, plane acoustic wave in the limit M -» 0. However, the author's previous 
asymptotic studies have shown that it is oblique acoustic waves which are primarily 
responsible for leading-edge receptivity at low Mach numbers. The approaches of 
Murdock and Reed are not easily modified to treat the case of obliquely incident 
acoustic waves. In addition, even for the parallel plane wave case, the actual unsteady 
flow far from the body consists of the superposition of an incident, plane acoustic 
wave and a scattered, outgoing, cylindrical acoustic wave. Since the amplitude of the 
scattered wave decays only as the square root of distance from the body nose, the 
assumption of a uniform axial pulsation at a finite distance from the body is question- 
able. 

Our computational approach takes advantage of singular perturbation concepts 
and has several advantages compared to the approaches utilized by Murdock and 
Reed. In contrast to these authors, we utilize analytical solutions for the steady and 
unsteady components of the potential flow outside the boundary layer. The singular 
perturbation approach allows us to correctly account for the scattering of the incident 
acoustic wave by the body surface, and to treat the important case of obliquely 
incident acoustic waves. Furthermore, with our singular perturbation approach these 
features can be captured without considering flow compressibility in the Navier- 
Stokes computation. This is discussed in more detail in the following section. Finally, 
in contrast to the approaches of Murdock and Reed, our computational domain can be 
restricted to the viscous portion of the flow field. This significantly decreases the 
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computational effort and provides advantages relative to the accuracy and resolution 
of the results. 

Two additional modifications which we have introduced into the Navier-Stokes 
code also merit discussion. The code is a velocity-vorticity formulation. The convec- 
tion and diffusion of the vorticity field is calculated in the first half of the time step. 
In the second half of the time step, an iterative solver is applied to the divergence and 
curl relations in order to update the velocity field. In our early unsteady flow calcu- 
lations, we experienced difficulties due to a gradual buildup of numerical errors 
related to incomplete relaxation in the iterative solver. This difficulty could be sup- 
pressed by significantly reducing the time step, but this resulted in prohibitive com- 
putation times. 

After some study, we realized that this difficulty was due to the way in which 
time variations in the potential flow field were being introduced into the code. The 
change in the inviscid slip velocity was being introduced through the outer boundary 
condition in the second half of the time step, and in each iteration of the velocity 
solver the adjustment in the potential flow field "diffused" inward only one grid level. 
Thus, a large number of iterations were required to "numerically diffuse" the change 
in the potential flow field. This characteristic of the numerical algorithm was incon- 
sistent with the inherent physics of the problem, since adjustments to a potential flow 
are instantaneously transmitted throughout the field via the unsteady pressure distri- 
bution. 

To more accurately model the physics of the problem, we modified the numer- 
ical algorithm such that the change in the potential flow field was introduced 
throughout the computational domain prior to applying the iterative solver in the 
second half of the time step. This modification dramatically improved the behavior 
of the algorithm. Much larger time steps can now be used with no indications of 


numerical difficulties. 
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After incorporating this change in the code, calculations were carried out over 
time intervals corresponding to several periods of the unsteady motion. The results 
obtained were in general very encouraging, but contained one feature which we view 
with some skepticism. Instability waves appear to be originating near the downstream 
boundary of the computational domain. It is possible that these instability waves are 
being generated by Morkovin's mechanism of unsteady free-stream pressure gradients, 
which could arise due to the interaction of free-stream acoustic waves with the curva- 
ture of the surface. To investigate this possibility, we are currently performing calcu- 
lations for the flat plate case, in which Morkovin's mechanism would not be operative. 
However, in recent work 15 we have shown that the instability waves produced by 
Morkovin's mechanism have exponentially small amplitudes whenever the length scale 
of the free-stream pressure gradient is long compared to the instability wavelength. 
Hence, it seems likely that another mechanism is producing these waves. In our 
previous computations, we found that the solutions are fairly sensitive to the down- 
stream boundary condition. Thus, we are also exploring the possibility that these 
waves could be spuriously generated by the downstream boundary condition. 

The asymptotic theory for leading-edge receptivity shows that the TS waves are 
generated near the leading edge where ux/U^ - 0(1). Hence, in a transient problem 
in which the unsteady motion is initiated at t - 0, TS wave motion should appear first 
in the upstream region. This TS wave motion should then gradually move down- 
stream as the unsteady flow field evolves. Thus, in the early stages of the unsteady 
flow development, the unsteady motion at the downstream boundary should simply 
consist of the Stokes wave induced by the local unsteady potential flow. 

In the computations described above, the time-harmonic motion sin ut was 
"switched on" at t = 0. However, we utilized the "steady-state" time-harmonic Stokes 
wave solution in the downstream boundary condition. The "switch on" of the time- 
harmonic Stokes problem has been investigated by Panton 1 *, who found that the tran- 
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sient died off in less than one period of the motion. However, we are concerned that 
the initial inconsistency in the outflow boundary condition may produce a spurious 
"numerical source" of receptivity. To investigate this issue, we are presently incorpo- 
rating the transient Stokes solution into our downstream boundary condition. 

The computational work described in this section is being closely coordinated 
with the development of the asymptotic theory, which is described in the following 
section. 


Asymptotic Theory 

We now discuss our progress on the asymptotic theory for leading-edge recep- 
tivity. The asymptotic analysis requires a sophisticated understanding of applied 
mathematics and fluid dynamics. Thus, in order to make timely progress on this 
project, a Post Doctoral student is required. Paul Hammerton was selected from a 
number of candidates for this Post Doctoral position. Hammerton received his Ph.D. 
in 1990 from the University of Cambridge. His dissertation advisor was David 
Crighton in the Department of Applied Mathematics and Theoretical Physics. 
Hammerton accepted a two-year appointment at the University of Arizona, and 
arrived at Arizona in early November. He is taking primary responsibility for 
carrying out the detailed asymptotic analysis of the viscous motion, and has rapidly 
come up to speed on the subject of boundary-layer receptivity. 

Both the asymptotic theory and the Navier-Stokes computations rely on analyti- 
cal solutions for the steady and unsteady components of the potential flow field. 
Thus, this part of the analysis was developed by the author prior to Hammerton's 
arrival in Arizona. The steady component of the potential flow field corresponds to 
inviscid flow past a parabola and is well known. The unsteady component of the 
potential flow field corresponds to scattering of a plane acoustic wave by a semi- 
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infinite parabola in a low Mach number flow. We have analyzed this problem using 
singular perturbation techniques. 

It might at first appear that the acoustic scattering problem in the limit M — * 0 
could be analyzed simply by assuming that the flow is incompressible. However, the 
author and Heinrich® showed that acoustic scattering problems involving semi- infinite 
bodies are singular in the limit M -* 0. Thus, the local flow field in the vicinity of 
the leading edge cannot be calculated without also considering a global region in 
which the acoustic wavelength A is the appropriate length scale. For the semi- infinite 
parabola with S = ujr n jU (X) = 0(1) and M = U^/c « 1, the parabola appears to be 
very thin when viewed on the scale of the acoustic wavelength A = 2ttc/w. Thus, the 
leading order term of the global solution is the classical problem of diffraction by a 
semi-infinite, zero- thickness plate. We utilize the Wiener-Hopf technique 17 to obtain 
the leading term of the global solution in the form 


4> 0 (x,y) 


Eac_i_i Jk (* cos » ac - y sin 9 ac ) 

pc \ k 


oo 

*<»c sgn(y) exp(->/A 2 - k 1 |y| - /A x) 

2 "v^J-oo " 


(A + k cos 6 ac ) \/A + k 


rfA je-' 1 ** 


(2a) 


Here (x,y) are dimensional coordinates centered at the focus of the parabola, p ac /pc is 
the magnitude of the velocity fluctuation associated with the incident acoustic wave, 
9 ac is the incidence angle of the acoustic wave measured with respect to the body axis, 
and k - w/c is the acoustic wavenumber. For matching with the local region near the 
leading edge, the asymptotic expansion of this solution as r - Vx 1 + y 1 -» 0 is 
required. From Fourier transform theory, it is known that the behavior of the 
integral in (2a) for small (x,y) is related to the behavior of the transform function for 
large A. Thus, we find that the small argument expansion of $ 0 (x,y) takes the form 
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<t> 0 (x,y) 


Pac _ I 

pc L 


. 3/2 e ,> / 4 ■ 9 a c r o 

2 ^ 7 = sin Vr cos x 

vir/c ^ 


+ cos 0 ac r cos 0 + 


0(r 3 /*) 


,-IWf 


(26) 


where ( r,6 ) is the polar form of ( x,y ). 

In the local leading-edge region, the appropriate length scale is the nose radius 
r n . The asymptotic analysis of the viscous motion utilizes coordinates nondimension- 
alized on the length scale U^/w, with the nondimensional nose radius S = VTn/U* 
assumed to be <9(1). Hence, anticipating the viscous analysis, it is convenient to also 
express the potential flow solution in the local leading-edge region in terms of coordi- 
nates nondimensionalized by Since r n « A, the unsteady component of the 

potential flow in the local leading-edge region is incompressible and consists of a sum 
of eigensolutions 

%{x,y) = + BRe[Vz - / Vs/5\ 

+ CRe[[\/Z- / Vs72] 2 ] + (3a) 

where Z = w(x + iy)/U ao and A, B and C are undetermined coefficients. Expanding 
(3a) for large Z and matching with (26), the constants are determined as 


A = 



B 


2 3/2 e<>/4 sin fl ac /2 

Vt Vm 


and 


C = cos 6 ac . 


(36) 


Physically, the function whose coefficient is B corresponds to an asymmetric, 
unsteady potential flow around the nose of the parabola, while the function whose 
coefficient is C corresponds to a symmetric, unsteady potential flow. For the case of 
an incident plane wave parallel to the axis of the parabola ( 0 ac = 0°), the coefficient B 
vanishes and C is 0(1). This is the case modeled in the computations of Murdock and 
Reed. In contrast, for the case of an obliquely incident acoustic wave, B = 0(M' X ^). 
Thus, at low Mach numbers, an oblique acoustic wave generates a strong asymmetric 
flow around the leading edge. This asymmetric flow in turn produces a strong 
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unsteady response in the boundary layer, leading to high values of leading-edge 
receptivity. Equations (3) provide the outer boundary conditions for both our 
Navier-Stokes computations and our asymptotic analysis of leading-edge receptivity. 

We now turn to the asymptotic analysis of the viscous motion in the boundary 
layer. As discussed above, the analysis of the viscous motion involves two distinct 
streamwise regions, the LUBLE region near the leading edge in which ux/U^ = 0(1), 
and the OSE region farther downstream in which uix/U^ = 0{e~ 2 ), where e = 
(uju/U^ 1 ) 1 ^ 6 is the small parameter in the asymptotic analysis. The asymptotic struc- 
ture of the unsteady viscous motion is illustrated in Figure 1. 

The viscous analysis is most conveniently developed in a parabolic coordinates, 
defined by 

* = TT 5 (** " ’J 2 ) ^ y * > (4a, b) 

where x and y are Cartesian coordinates in the streamwise and cross-stream direc- 
tions, respectively. We have nondimensionalized the physical coordinates (x,y) by 
Oqo/w, so that the LUBLE region corresponds to ( = 0(1) and the body surface is 
given by f) = S 1 ^ 1 . To thoroughly investigate the effect of nose bluntness, the analysis 
is developed for the case S = 0(1). The viscous motion in the boundary layer is 
driven by the unsteady slip velocity associated with (3). Expressing (3) in parabolic 
coordinates and setting 1) = S 1 ^ 1 , the unsteady slip velocity is given by 


Pac [ 2 e**/* 


It can be seen from this expression that the nose radius of the parabola limits the peak 
velocity associated with the strong asymmetric flow due to obliquely incident acoustic 


waves. 
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Since the thickness of the mean boundary layer is of the order of the square 
root of the Reynolds number, we introduce a boundary layer coordinate rj where 

T) = S l/1 + e 3 rj. (6a, b) 

Separating the streamfunction into a steady flow component $(£,»?) and an unsteady 
disturbance rp((,T])e ' iM , and linearizing the equation for 0, the unsteady flow in the 
boundary layer is found to satisfy the linearized unsteady boundary layer equation 
(LUBLE) 


^ r lV T t + ^>TV + ^5 ~ ^r) ~ $1/ 


+ 



S t d_ 

~h* ~ 7r d$ 


u sl(0 


(7) 


where H 2 = ( 2 + S arises from the Jacobian of the transformation from Cartesian to 
parabolic coordinates. The right hand side of this equation is the unsteady pressure 
gradient impressed on the boundary layer by the unsteady potential flow. At rj = 0, 0 
and rp , , are zero, and rj> n must match with the inviscid slip velocity as rj — ► oo. 

For £ = 0(1), solutions for and 0 can be found only by using numerical 
methods, which we shall discuss in due course. However, for large values of £ well 
downstream of the nose, an asymptotic solution for (7) can be derived. This asymp- 
totic solution consists of a Stokes wave which is related to the inviscid forcing field, 
plus an infinite set of asymptotic eigensolutions. The first asymptotic eigensolution is 
the precursor of the TS wave, and thus is of fundamental importance in the theory of 
leading edge receptivity. 

In calculating the asymptotic eigensolutions of (7), we utilize Van Dyke's 18 
result for the mean boundary layer far downstream on the surface of a parabola in a 
uniform stream. He found that the streamfunction for the mean flow has the expan- 


sion 
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* = £ [ Fin) + ^ log (p/S) a, [ nr - F] + ±[b 1 [nF - F] + /(»,)] + ... J (8) 

where A r and B 1 are numerical constants, F is the Blasius function and / satisfies a 
prescribed ordinary differential equation. This expansion shows that, far downstream, 
the boundary layer on a parabola approaches that for a flat plate. Physically, the term 
involving B x represents a slight shift in the effective origin for this flat plate. 

Utilizing (8) for the mean flow description, we develop an asymptotic solution 
of the LUBLE (7) for { » I. The most important feature of this solution is the first 
asymptotic eigensolution. The asymptotic eigensolution has a two- layer structure in 
the direction normal to the surface. The inner layer corresponds to 17 = while 

the the outer layer has the same thickness as the mean boundary layer, t) = 0(1). In 
analyzing the inner layer, it proves convenient to introduce a modified transverse 
coordinate m = £r;(l + S/ 2£ 2 ). The use of this coordinate simplifies the solution by 
taking partial account of the parabolic geometry. In the inner layer, the expression 
for the first asymptotic eigensolution has the form 


= <^(5) £ 2t [p 0 (m) + f 3 Px ( m ) + •••] (9 a) 


where 


no = 


3 U' 0 

. 


1 + 3^ ^ In (S/0) + 3 



2A, 



(9b) 


and {7J, denotes £/'( 0). Equation (9) generalizes Lam and Row's 4 asymptotic eigensolu- 
tion for the flat plate to the case of a parabolic cylinder. The constant A in (96) is a 
(real) eigenvalue entering the governing equation for p 0 (m). The leading edge blunt- 
ness has produced logarithmic terms in the exponent T(0- The value of r, a numeri- 
cal constant, is fixed by the solvability condition for pjm), the second term in (9a). 
This gives 
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T = T s = o + 


S X 


2 u: 


7T 


i [C 

^ j°° m w 2 dm 


(10a) 


where r 5=0 is the value for the flat plate case and w(m) = p„(m) is expressible in 
terms of an Airy function. Thus, in the exponent r, the leading edge bluntness 
produces an order-one change to the flat plate result. Goldstein 6 calculated the flat- 
plate value r 5 _ 0 ~ -0.69. Evaluating (10a), r has the complex value 


t at - 0.69 - S 0.51 (1 - 0 


(106) 


for the case of a parabolic leading edge. The imaginary part of r is allied to the 
appearance of logarithmic terms in T(£). Physically, the imaginary component of r is 
equivalent to a logarithmic term in the phase of the asymptotic eigensolution. 

The coefficient C : multiplying the asymptotic eigensolution in (9) is arbitrary 
within the asymptotic analysis, and hence must be found from a numerical solution of 
the LUBLE. This is discussed below. In general, C t is a function of the nondimen- 
sional nose radius S and of the specific free-stream disturbance. 

The coordinate transformation (4a) shows that distance downstream is propor- 
tional to £ 2 . Hence, the cubic term in (96) implies that the wavelength of the asymp- 
totic eigensolution is proportional to jc" 1 / 2 . As the wavelength decreases with distance 
downstream, the pressure fluctuations induced by the unsteady displacement thickness 
of the asymptotic eigensolution increase in importance. When £ = this self- 

induced pressure gradient exerts an order-one influence on the unsteady motion. 
Since the two-layer solution (9) ignores this pressure field, it becomes invalid in this 
region. Therefore, the two- layer solution must be replaced by a new asymptotic 
solution. This new solution turns out to be the high Reynolds number, small wave- 
number asymptotic expansion for the Orr-Sommerfeld equation (OSE). This asymp- 
totic expansion corresponds to a triple-deck structure, in which the two layers of the 
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previous solution are supplemented by a third layer just outside the mean boundary 
layer. This third layer is the potential flow which enters the interactive pressure- 
displacement relation. 

Anticipating the asymptotic matching of the LUBLE and OSE regions, the TS 
wave in the OSE region has the expansion 



r 

r ?i l 

\p = r 2r 7 (»j,e) exp 

i 

€ 3 ^ 

di\ 


where ^ = «£ = 0(1), t is given by (10), is a slowly varying amplitude, -y is the 
mode shape and k is the complex wavenumber. With suitable changes in the defini- 
tion of the wavenumber a, phase speed c and Reynolds number R, the OSE for the 
parabolic geometry can be reduced to the standard form for the flat plate. Specifi- 
cally, we set 


Q - Q s=o 


,(S f 


1 c - C S=0 


1 + 5 


2*i 


and 


aR = (aR) s _ 0 


l+ S W 


(12a,6,c) 


Here the subscript 5=0 refers to the appropriate values for the flat plate case. 

At first sight these corrections would appear to have little effect on the form of 
the TS wave structure, beyond modifying the slow spatial variation of the wave amp- 
litude. However, a careful analysis of the modified equation demonstrates that the 
standard results of the triple-deck theory for the flat plate case are not directly appli- 
cable. The mean flow pressure gradient leads to a crucial difference in the boundary 
layer profile close to the body surface, namely t/"( 0) * 0. Thus the complex wave- 
number obtained by Goldstein 1 for TS waves on a flat plate can not be readily 
modified to the parabolic geometry. However using similar methods, we have 
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obtained the following expansion of the wavenumber for the parabolic cylinder. 


*(£ 1 ) = [k 0 + £ Kj + € 2 /c 2 ] + 5 [e 2 In e k L<2 + < 2 ftj + Q(e 3 In e) , (13a) 


with the e 2 In e term arising due to the first perturbation term in (96). The functions 
k 0 , /Cj and /c 2 are the flat plate values given by Goldstein 2 . The first two terms gener- 
ated by the parabolic geometry are given by 
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H(z) = gl ' Sy/2 diM . 

J1 
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r/T(y) 


(136) 


(13c) 


(13c/) 


In these expressions, Tis defined by H(Y) = ^ where ^ = iJU' Q and A 1 and B 1 are the 
numerical constants appearing in (8). It is interesting to note that the G(e 2 In e) cor- 
rection due to the parabolic geometry precedes the first correction due to nonparallel 
flow effects which is 0(e 3 ) in Smith's 3 flat plate theory. 

The higher order terms for k have also been determined, up to the e 3 term 
which also determines the slowly varying amplitude function A(^ 1 ). This solution in 
the triple-layer region has been rigorously matched back to the asymptotic expansion 
(9a) for the LUBLE. This matching shows that the coefficient in (9a) determines 
the amplitude of the TS wave, and hence is the "Receptivity Coefficient." The further 
development of the TS wave with distance downstream over the parabolic cylinder, as 
described in (13), is not central to understanding the actual receptivity process. 
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However, the detailed effect of the parabolic geometry on the TS wave must be fully 
determined in order to allow comparison with the ongoing numerical investigation of 
Dr. Gatski. 

As an illustration of these effects. Figures 2 and 3 show preliminary compari- 
sons between the TS wavenumbers for a parabolic cylinder with S = 1 and for a flat 
plate, with e = 0.1 in both cases. The normalized coordinate ^ has absorbed some of 
the effects of the parabolic geometry, and hence the neutral stability points for the 
two cases in Figure 3 are not markedly different. However, the favorable pressure 
gradient over the parabolic cylinder modifies both the wavelength and the growth rate 
of the instability wave. The favorable pressure gradient decreases Re(/c), correspond- 
ing to an increase in the TS wavelength (Figure 2), and further stabilizes the boundary 
layer, as is evident from the growth rates plotted in Figure 3. 

At small £j, the asymptotic results plotted in these figures become invalid as the 
region governed by the LUBLE (7) is approached and geometric effects become 
dominant. The lower limit = 0.1 corresponds to £ ~ 1.6 which is well within the 
region of significant wall curvature. 

To determine the Receptivity Coefficient C lt a numerical solution of the 
LUBLE for the particular free-stream disturbance of interest is required. Essentially, 
the numerical solution of the LUBLE for £ » 1 is compared with the expression (9a) 
for the first asymptotic eigensolution. We are now beginning to address this task, and 
are considering the merits of various numerical schemes. An important aspect of our 
approach is to utilize our knowledge of the analytical structure of the large £ behavior. 
A critical feature of the problem is that for real u and £, the first asymptotic eigenso- 
lution makes only an exponentially small contribution at large £. Thus, for real £, it 
would be extremely difficult to determine C x from a direct comparison of the numer- 
ical and asymptotic results. However, the LUBLE is analytic in £, and hence we 
extend the problem into the complex £ plane and calculate the numerical solution 
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along a ray where the first asymptotic eigensolution is dominant. As in previous work 
by the author, we shall analytically factor out the now rapidly growing exponential 
factor e' r (£) of the first asymptotic eigensolution before discretizing the problem. 

We anticipate that the development of a computer program to solve the LUBLE 
and extract the receptivity coefficient will take only a few months. By the time of 
the ICASE/LaRC workshop this summer, we should be well along in a parametric 
study based on our asymptotic theory. During the workshop we plan to make exten- 
sive comparisons of our asymptotic results with the numerical results generated by Dr. 
Gatski's Navier-Stokes code. 


Concluding Remarks 

In the first year of this research program, significant progress has been made on 
both the asymptotic theory and the collaborative Navier-Stokes computations of 
leading-edge receptivity. This project complements the NASA LaRC in-house 
activity on boundary-layer receptivity, which is focused mainly on localized rather 
than leading-edge mechanisms and involves the author's ex-student. Dr. M. Choud- 
hari. We have developed a strong, mutually beneficial collaboration with Dr. Gatski 
in Navier-Stokes computations, and foresee additional opportunities for collaboration 
with other NASA computational fluid dynamicists. In particular, the author and Dr. 
C. Streett recently discussed possible interactions on spectral calculations of leading- 
edge receptivity, and we plan to explore this further during the ICASE/LaRC 
workshop this summer. Additional opportunities for collaborative activity would arise 
in the second phase of this research program, which extends the investigation to 
supersonic Mach numbers. 

This project involves subtle physical mechanisms and sophisticated mathemati- 
cal techniques, as can be seen from the above discussion of our progress on the 

, t! 
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asymptotic theory. The significant background work in this area has been performed 
by outstanding senior researchers such as M.E. Goldstein and F.T. Smith. Thus, when 
the author proposed this project, he emphasized the importance of utilizing a Post 
Doctoral student. It would be quite difficult for a beginning graduate student to work 
through the relevant background literature, and even harder for him to understand 
how these analyses could be significantly extended. Indeed, if this project was 
attempted with a beginning Ph.D. student, the author would anticipate very little 
progress for at least a couple of years. The low Mach number study alone would be 
appropriate for a three (or four!) year Ph.D. project. 

The argument is sometimes advanced that projects involving graduate students 
are preferable due to their educational component. Post Doctoral studies also have a 
strong educational component and in fact are becoming increasingly important for in- 
dividuals seeking an academic position at a first-rate university. The education 
provided by Post Doctoral study is equally valuable for future positions in government 
or private research laboratories. An alternative viewpoint is that the choice of a 
Master's, Ph.D. or Post Doctoral student should be based on the difficulty of the 
project. In all three cases, the project contributes to the education of the student and 
to his future career development. 

In summary, at this time our research on leading-edge receptivity is definitely 
on schedule, and possibly ahead of schedule. The asymptotic theory for the case of a 
symmetric, low Mach number mean flow is nearly complete. The only remaining task 
is the development of a computer program to solve the LUBLE and extract the recep- 
tivity coefficient C r The computer programming should take only a couple of 
months. During the ICASE/LaRC workshop this summer, we shall utilize the asymp- 
totic theory to perform a parametric study of the influence of nose bluntness on 
leading-edge receptivity for low Mach number flows. Extensive comparisons of the 
asymptotic results and Dr. Gatski's Navier-Stokes simulations will also be carried out 


during the workshop. 
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FIGURE 1. The asymptotic structure for the unsteady disturbances in the boundary layer 
of a parabolic cylinder. 
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FIGURE 2. Variation of the real part of the TS wavenumber (multiplied by e) as a function 
of the downstream distance (Re(«) vs. = (fi/££') 2 ) for e = (wi//f7oo 2 ) - ? = 0.1. The 
flat plate result is compared to that for a parabolic surface with s = ur n /U<x> = 1.0, where 
r n is the nose radius. 




Growth Rate 
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FIGURE 3. Variation of the growth rate of the TS wave (multiplied by e) as a function 
of the downstream distance (— Im(/c) vs. $ = (ti/U^) 2 ). Flat plate results are compared 
with those for a parabolic cylinder with s = 1. The stabilising effect of the favourable 
pressure gradient can be seen, this effect decreasing with distance downstream as the body 
curvature approaches zero. 



